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Abstract 

The statistical characteristics of wind-shear turbulence are modelled. 
Isotropic turbulence serves as the basis of comparison for the anisotropic 
turbulence which exists in a wind shear. The question of how turbulence 
"scales" in a wind shear is addressed from the perspective of power spectral 

density. 
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Nomenclature 


of poor £'jauty 


A (r,t): 

U - 

C (r,t): 

l) - 


r : 

k 

r r — ) f r (*••): 

1 2 

J ( r\ t ) : 

ij - 

M ,M : 

1 2 

P (q, r ) : 

m 

Qlq. r, t ) : 


R(r.t): 


S ( r , t ) : 

1 jm ~ 

U : 

l 

0 : 

\ 

V: 


a, b: 


flr\t)/g(r,t): 


k ( r , t ) : 


r : 

r : 

1 1 : 
l 

w ( x , t ) : 
x , z : 

A ( I ) : 


second-order anisotropic tensor 

two-point time-dependent second-order velocity correlation 
tensor 

decay constant defined in Equ. (12) 
polynomial basis functions for Q( • “ * ) 
second-order isotropic tensor 
mean-flow wind shear gradients 

homogeneous polynomial of order m in the variables, q»r 
third-order three-point scalar correlation function 
normalized scalar velocity correlation function 
third-order two-point velocity correlation tensor 
total flow velocity 
mean-flow velocity 
airplane flight speed 

time-dependent parameters characterizing anisotropic 
turbulence 

longitudinal/transverse correlation functions for isotropic 
turbulence 

scalar function peculiar to third-order velocity moments of 

turbulence 

separation vector 

magnitude of separation vector 

turbulent flow velocity 

turbulence downwash 

Cartesian coordinate variables 

integral scale for isotropic turbulence 
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AA( t): 

A(fc, t): 
a: 

ft, H ' 

K : 

i : 

<r( t ) : 

A<r( t ) : 
w: 

0 : 

0( K, t ) : 

0 ( K , t ) : 

0 U) ( K , t ) : 


change in integral scale due to anisotropy 

energy-scaling factor 

airplane glide-slope angle 

dimensionless constants 

non-di mens ional ized wave-number 

imaginary element 

isotropic turbulence intensity 

change in intensity due to anisotropy 

angular frequency 

dummy variable 

power-spectral density 

power-spectral density for isotropic turbulence 
anisotropic part of power-spectral density 
time-lag 
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1 . Background 

Meteorological events such as thunderstorms and unstable frontal systems have 
long been considered dangerous from the perspective of aviation safety. 
Investigations into the alarming number of recent low-level aircraft accidents 
involving thunderstorms around airports have revealed one meteorological 

culprit a strong blast of air directed toward the ground which has come to 

be known as a "microburst". 

The mechanism by which the microburst is created is a complex and sometimes 
violent one. From the investigation of a particular airliner crash at JFK 
airport, much was learned about the convective atmospheric dynamics which 
cause the microburst. It was postulated by Fujita* that a violent atmospheric 
disturbance may occur when a moist upper air fract drops precipitation through 
a relatively dry layer below it. As the precipitation evaporates, the dry 
layer cools; consequently, a stream of dense cold air suddenly replaces the 
ground-heated low-altitude air. Furthermore, as the resulting downdraft 
penetrates through the existing horizontal flow, turbulence is invariably 
generated. The downdraft eventually impacts the ground and spreads out 

radially. During the evolution of a microburst, the dimensions of the 
low-altitude phase may be on the order of 15-30 meters high with a radial 
spread of 2.5 kilometers. Realization of its relatively small meteorological 
size serves as the motivation for the term "microburst". Wind velocity 

information subsequently obtained from flight recorders has sustained Fu.jitn's 
explanation. The typically distinguishing characteristics of the microburst 
are winds which are localized and variable. Abrupt changes in wind speed 

and/or direction over a short atmospheric distance are referred to in the 
literature as "wind shear", a term coined by Fujita himself. For purposes of 
this study, however, wind shear will be defined as a flow which lias a 
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spatially nonconstant mean-velocity profile. Flight through a "ulnd 

shear-infected area" will tax both pilot and aircraft since related velocity 
gradients may drop airspeed to critically low levels. In summary, microbursts 
and wind shear pose an immediate threat to flight operations at low altitudes 
due to their adverse effects on aerodynamic lift/aircraft response. 

The scenario of microburst encounter by an aircraft in the landing mode is 
shown in Figure 1. First the aircraft sense an uncharacteristically strong 
head wind. Although turbulence typically accompanies the increased airspeed, 
pilot confidence is high during this phase of the encounter since aircraft 
performance increases with the additional aerodynamic lift. Because the 
intent is to land the aircraft, the pilot intuitively trims. Eventually, the 
head wind becomes a tail wind, and at this point performance is immediately 
and seriously degraded since the distinct head-to-tai l wind swing 
airspeed enough to possibly lose flying speed. In the last ten years, 
estimated 575 people have been killed in commercial aircraft accidents in the 
United States alone due to the microburst phenomenon. Likewise, 

■inroburst-induced accidents have also been reported in Great Britain, 
Germany, France, Italy, Australia, and Japan. Consequently, a fundamental 
understanding of microbursts, related wind shear and turbulence, is a 
h 5 gh-priori ty research issue. 

Flight science authorities agree that the "solution" to the wind shear 
"problem" is multi-faceted. The existing research programs of NASA and FAA 
focus primarily on three elements, II hazard characterization, 21 sensor 
technology, and 3) flight management and operations. Hazard characterization 
is the study of the physics of the microburst phenomenon. Inherent within 
this phase of the research effort is analysis of aircraft aerodynamics in wind 
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Wind shear velocity profiles, rainfall effects, and 


shears and in heavy rain, 
turbulence models are important contributions from this component of the 
programs. The sensor technology component deals with the prediction and 
detection of potentially dangerous meteorological events. This component 
subdivides into airport-fixed and airborne sensor technologies. Next 

Generation Radar fields ( NEXRAD ) are presently being investigated by NOAA, FAA 
and USAF for future use as airport radars. The other detection device is the 
in situ look-ahead sensor. This onboard computer system is being devised to 
scan the forward atmospheric environment, evaluate a "hazard index", and 
annunciate a warning to the pilot in a time period adequate to ensure either 
avoidance or escape of the threat. Look-ahead sensor technology is beyond 
state-of-the art, and is not expected to be operational until the mid-1990’s. 
The third component of microburst research is flight management. Simply put, 
this phase concentrates on informing the pilot on how to safely get the 
aircraft out of a microburst encounter subject lo aircraft performance 
constraints. An important element of such studies is flight simulation. 
Simulators allow pilots to experience threatening wind shears in a controlled 
environment, with intent to better prepare them for potential real-life 
events. Realistic wind shear representations in flight simulators are 
practical and economical, and can help crews coordinate their escape efforts 
in critical situations. Realism in the simulation of flight through hazardous 
atmospheric environments has improved with the introduction of lateral and 
vertical winds, as well as vortex and turbulence influences. The focus of the 
present work is in the hazard characterization aspect of the programs: its 

principle objectives are to mode) the fluid-dynamical characteristics of 


turbulence associated with low-level sheared mean-flow, and quantify the 
effect of these characteristics on aircraft response. All of these issues 
will be herein addressed from the perspective of power spec tra d, dens i ty (psd). 
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2. Introduction 

Since microbursts pose a hazardous threat to flight operations only at low 
altitudes, the role that the ground plays as a fixed boundary is therefore 
important in terms of turbulence modelling. The presence of the ground 
generates in the microburst a brand of turbulence which when sensed by the 
aircraft, during its take-off or landing is both nonstationar y and anisotropic - 
The objectives of this investigation are thus to model the statistical 
characteristics of same, and gain a quantitative understanding, in terms of 
scaling effects, of the stochastic nature of the turbulence sensed by an 
aircraft during a microburst encounter. 


The task of mathematically describing the turbulence associated with a typical 
microburst (shown in Figure 1) is straightforward. The pioneering works of 
Taylor 2 . Kerman and Howarth 3 , and Robertson 4 have illustrated the well-defined 
fashion in which the concepts of two-point velocity correlations can be used 
to quantify the kinematical aspects of isotropic turbulence; Batchelor’ 
extended these concepts to the case of ax i symm etric, turbulence that which 
has a defined form of anisotropy. Axisymmetric turbulence, with its tendency 
to prefer one spatial direction, is the logical candidate for modelling wind 
shea,- turbulence (with its particular anisotropy resulting from the microburst 
impacting the ground). Although classical axisymmetric theory produces 
turbulence velocity correlation functions that are posi t ion- i ndependenl 
(homogeneous) and time-varying, the turbulence under consideration is assumed 
inhomogeneous and in a steady-state ("frozen") condition: nonetheless, it can 
be easily shown 6 that correlation functions for frozen. inhomogeneous 
turbulence may be converted into time-evolution equations for homogeneous 
turbulence by duly scaling out a characteristic velocity. Accordingly. 
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r (x.r) - C ( t , r ) -» C (r, t) 
i j -- • * ) " 1 J 

physics of the inhomogeneous, 
of classical turbulence; the 
the accepted standard. 


and S (x.r) ~ S it, r) S 

tjm - ~ ijm 

steady-state turbulence aligns 
(r, t ) -notation rather than the 


( r , t ) , and 1 he 

imj ~ 

with the physics 
( t , r ) -notat ion is 


Recent work 7 , which incorporates axisymmetric theory, has resulted in the 
two-point velocity autocorrelation for anisotropic turbulence downwash. viz. 


£ f r t) = < w ( x . t. ) w ( x + r\t)> - {cr^g} + {a|r| + b} ; 

33 


( 1 ) 


here a ~ 0(a). where a is the airplane glide-slope angle, b is a parameter to 
be defined later, and g is the transverse correlation function of isotropic 
turbulence. Equation (1) is based on the assumption that during take-off or 
landing the aircraft senses only a small amount of anisotropy in the 
(necessarily) large amount present in the microburst. Note that the velocity 

correlation separates into isotropic and anisotropic components, and takes the 

fisotropic \ + /anisotropic! where 
general form, c 33 ( r ’ tJ ~ /correlation/ /correlation) 


fisotropic 
\ cor re lat i 


on 


2 

= cr g 


and 


{ anisotropic 
correlat 


pic l = a 
ionj 


r I + b. 


The contributions of isotropy and anisotropy to the total correlation are 
shown in Figure 2. 


An important aspect of the ensuing analysis 
microburst. A common correlation used for 
uniformly modulated, viz. fdr.t) — (t)R(r) 

modulator; implicit, in this correlation is 


is how turbulence "scales" in a 
nonstationary processes is the 
where .r 2 ( t ) is the time-varying 
1 he constraint imposed on the 


integral scale — that of 
general form has been 
intensity and integral 


time- invariance. In recent studies', however, a more 
formulated which incorporates time-dependence in both 
scale, viz. C(r.L) = tjR(?D where v = T»(t) = <r°( t > an,J 
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g = f;(t) = 1 / A ( t ) . Allowing the correlation to evolve in such a manner 
permits its geometric shape to distort in time, while its r-algebraic form is 
retained. Figure 3 illustrates the "self-preserving" behavior of such a 
correlation function. With C(O.t) s <r = (t). intensity is clearly an 
ordinate-axis scaling effect while the normalized correlation function, Pdf;r>, 
automatically incorporates the integral scale/abscissa-axis scaling effect. 
Although self-preservation is introduced here strictly from the standpoint ol 
"scaling-effects", its hypothetical existence and subsequent analysis can be 
be .justified on purely "dynamical system" grounds (see Appendix). 


3. Power Spectral Density Analysis 

The power spectral density is a common way of illustrating energy distribution 
versus frequency, and is an immensely popular tool in both gust analyses and 
turbulence simulations: mathematically, it is the Fourier transform of the 

autocorrelation function. The psd for the aforementioned downwash formulates 
as 

<p{ w/V, t ) = J <w(x,t)w(x + r, t )>e ' dr = X {tr g + a |r ) + b)e dr, (2) 


where i 2 = -1. Since the correlation function consists of isotropic and 

anisotropic components, so too will its transform: specifically <p U,t) = 

0 m U.t) + where 


( 1+k‘)‘ 


(a) ,2. TaaI 1 -cos (Kb/a A) 

</> V.t) = 4<r A {--V — • <•*» 

J i Kb/ a A ) 


ll er -e ,/,(«. t) is the total psd sensed by the aircraft as it flies through the 
microburst, k = wA/V a non-dimensional i zed wave number. <r the iKolropic 
intensity, A t lie longitudinal correlation length (integral scale), and AA (lie 
change in integral scale due to anisotropy. The scaling effects discussed 
earlier are evident upon inspection of Equ. (3). 
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The distinct advantage that the psd representation has over its transform 
partner, the autocorrelation function, is its more defined statement of the 
modelling parameters of anisotropy, viz. M. a. and b. in relation to total 
turbulent energy. For a glide-slope which is maintained during a microbursi 
encounter, it follows that a * constant since aitl - Q(«). Hence, b plays the 
more important roie in terms of modelling the anisotropy due to the boundary 
presence. Since b is related to intensity by e 2 * b = ( <r * ho-| 2 , 4<r - change 

in intensity due to anisotropy, for neglected higher-order terms b becomes 


b = <r 


* [*}} * - 1 ") - ^ 


(4) 


One way to quantify the collective effects of AA, a, and b on the sensed 
turbulence is to introduce an "energy scaling factor", defined as AU.t) 
0 , 8 , ( K .t)V ,, (K.t). To analyze the effect of AA/A on AU.t), set b/(aA) 

, , „f A = 1000 f t , V = 225 ft/s and a ~ 0(a) 

constant; typical values oi <r - 5 it/s, A uuu 

= 0(3°) = 0(0.052 rad) are chosen, and Atr/cr is allowed to vary up to 0.5 (a 
feasible value for wind-shear turbulence severity). Incorporating Equ. (4) 
gives b/(aA) - 0(0.5). Figure 4 shows that for AA/A = 0.1, the energy scaling 
factor is as high as 35% for 0 =s k * 2 (0 * « - 0 . 45 rad/sec) -a frequency 

range which falls within that specified for the guidance and control phases of 

an aircraft mission involving a wind shear. Hence, a realistic estimate of 
the energy transferred to the airplane requires that anisotropy be accounted 
for in turbulence models and related flight simulations both. 


4. Constraints Imposed by Navier Stokes 

Knowledge of the sheared mean-wlnd/turbulence Interplay retires that the 
governing equation of the fluid be considered. The Havler-Stokes equation for 
a turbulent , high Reynolds number, steady-state fluid with both body and 
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pressure forces neglected is 


d ( U U ) 

t m 

dx 


^ Q. 


(5) 


wuh h = (J + u . the two-point correlation method yields 
t i i 


Trr C + « — C 
dt 1J <3x mj 


au ' 

2 — S + — — C =0 
5r m iraJ dx ' ,m 


m 


where C 


i 1 


C (r, t ) = <u U ' > and S 

' l J i j * m 3 


S I r , t 

ImJ ~ 


( 6 ) 


<u u u ' >; index 
1 m J 


contraction results in 


au 


au 


— — C + C + ~r C - 2 — — - — S - 0. 

dt il dx ml ., ml dr 1 m i 

m dx m 


( 7 ) 


Equation (7) requires the solution of separate isotropic and anisotropic 
dynamical problems, viz. 


d 1-2 -4 — S . = 0 


( 8 ) 


dt u 


dr iml 


and 


dU dU 

JL_ A + C + -r C = 0 

dt il dx mi ~ ml 


(9) 


The basis for a scaling comparison (the isotropic 
case ! will be considered first. 


7.10 , , . _ r 

where I + A — U 

wnei e ■ t } i} i j 


Introducing the f and g functional representations for the correlation tensors 

Hives a result whose first integral is the K^rman-Howart h equation for high 
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Reynolds number isotropic turbulence , viz. 


dt 


( <r f ) 


4k f dk I 
r dr J 

V. J 


( 10 ) 


where k = k(r.t) is the scalar function required to specify third-order 
velocity product moments. Since turbulence self-interaction is not altered by 
self-preservation, it follows that k(r.t) also has such an evolution with 
time. i.e. k(r.t) - k(^r) = klr/A). Incorporation of this feature into the 


1 1 



'-integration of the KArm&n-Howarth equation results in 


-1- {<r 2 (t)A(t)f = o- 3 C 
dt K 


(in 


where 


( 12 ) 


r - r /iiSliLL + lliifl d0 = negative scale-independent constant, 

c * ; l © a© J 

and 0 = r/A: a specific numerical value of C r will he determined later. 

Solution of Equ. (11) gives the time-dependence. 


C . 

Alt.) = — - — f c r ( t ) dt , 


( 13) 


<r 2 ( t ) 


for the integral scale. The power 

t n , n = ±1, +2, ... , converts Equ. (13) into 

2C 


law form for isotropic intensity, cr"(t) 


A(t) 


3n + 


k ^ ( n + 2 ) /2 ^ _ + ^ +2^ . ; 


( 14) 


1 / 2 , 


for n = - l , A( t ) ~ 0( t ) 


To consider the anisotropic part of the problem, viz. Equ. (9), gradients must 
first, be defined. The two wind-shear profiles chosen for analysis are 
illustrated in Figures 5 and 6. and will hereinafter be referred to as "Shear 
Flow 1" and "Shear Flow II", respectively. Note that both profiles have a 
"head-to-tai I wind swing", a phenomenon frequently observed in microburst. 
encounters. The mean velocity for Shear-Flow l is U ) = (U(z), 0,0} = {M^.0,0} 
as shown in Figure 5; H ( is negative on the headwind side of the microburst 
and positive on the tailwind side. Analysis of Shear Flow II is restricted to 
the two-dimensional case with U = (U(x) , 0, W(z) } and gradients of the form 


dt) 


dW 


dx dz 

Equ. (9) yields 

d 

"TTF 


= M > 0. 


Incorporation of Shear-Flow I characteristics inti 


( a I r I +b ) - M <r 2 <xr -4^- + M arc- 


dr 


2 3(f + 2g ) 
Or 


- 0 


( 15; 
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integration in r from 0 to w produces 

— — {o- 2 (t)AA(t)> - M<r 2 (t)aA(t) = 0 , 
dt 1 

00 

where J (a|r| + b) dr = </'AA(t). The solution to Equ. (16) is 


AA( t ) = — '■ J <r"(t)A(t)dt. 

<r 2 (t) 

For <r 2 (t) ~ t" and the integral scale of Equ. (14), Equ. (17) becomes 


AA( t ) 


4M C a 

1 k 

(3n + 2)(3n + 4) 


(n+4)/2 


n = ± 1 , ±2 , 


( 16) 


( 17 ) 


( 18) 


specifically, AA(t) ~ 0(t 3/2 ) for cr 2 ( t ) ~ t“\ A similar analysis for 

Shear-Flow II yields 


AA(t) 


4M C 

2 k (n+4)/2 

(3n + 2 ) ( 3n + 4) 


n = ±1, ±2, 


Like Shear Flow t, Shear Flow II has AA(t) 


0(t 3/2 ) for <r Z (t) ~ t 


5. Determination of C 

k 

The theory of turbulence strictly prohibits the presence in the analysis of 
any adjustable constants; therefore a full closure of the dynamical equation 
is not achieved until a value for C is determined. The fact that C is 

k K 

formed from an integral of R( ) (cf. Equ. (12)) explicitly suggests that 

aside from some fundamental limiting constraints, viz. r > oo (kl***)) ••• 0(r) 
and r 4 S { k ( — )} = 0, the exact functional form of k ( - - - ) is not critical, 
i.e. k( - • • ) itself is not necessarily a unique function. All that is required 
for any and all possible k ( - - * • , say k ( * - - ) — * k^l • * • ) , N = 1,2.3 is 

CO 

the! (4 f r ' k tr.tldr) = C be independent of N A preliminary step is 
N k 

o 

therefore to develop a representative algebraic form for k(r.t) k(r A). 

To accomplish this while simultaneously saying as little as possible about 
)/(...) itself, it is here proposed to kinematically describe the three-point 
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corre I at ion, 


<u(x.t)u(x + q, t )u(x + r,t)> = Q(q,r.t), where u(x,t) is only one 
component, of the three-component turbulence velocity field and is 

one-dimensional position vector, i.e. x = (x.0.0). etc. Note that for 

self-preservation. this correlation accordingly scales as U'q.r.tl 
<r J (t.)s(q/A.r/A). where s(0,r/A) = k(r/A); furthermore, the functional form of 

Q( ) is such that it must satisfy certain symmetry properties, viz. Qfc|, r*. 1 1 

= Q(r.q,t.) = Q(-q, r-q, t,) 11 - Since Q( ■ - * ) expanded in a Taylor’s series about 

q = p = 0 produces 

Q( q, r, t ) * a Q + a ? q + hr + a^ 2 + c^r + b^ 2 + • ■ • . 

it follows from the required symmetry that the time-dependent parameters, a^ 

b , and c , m = 0.1,2,..., cannot be completely arbitrary. In fact they must 
m 171 

1 O 

each be such that eventually Q( * * - ) expresses as 

00 

Qfq.r.t) = E ^ 

^ M u m m 


m = 0 


where each P (••■) is a homogeneous polynomial in q. r and individually must 

m 

themselves satisfy the same symmetry as Q(q,r,t). i.e. fMq> r > = P m (r,q) 

P ( _ q ^ r — q ) t m = 0,1,2,... • Fortunately, there are only two (2) fundamental 

such forms, viz. F/q.D = q 2 - qr + r 2 and F^q.r) = 2q 3 - 3q 2 r - 3qr 2 + 2r 3 , 

since any other higher-order such form can be represented as a combination of 

F (•••) and/or F (•••): there is no linear homogeneous algebraic form in q,r 

\ " 2 

which satisfies the required symmetry. For example, the fourth or dci 
polynomial is identical to {F^q.rl} 2 , while the fifth-order one is identical 
to (F 1 q , r ) F tq,r)}. The existence of such a complete set of polynomial basis 
functions for any sequence of systematically generated algebraic forms was 
original ly established in a theorem formulated by David Hilbert m 1890 . 

The theorem need not be quoted here, but the essence of it is that if an 
infinite sequence of algebraic forms (the sequence h_re is P () t !• 

P ( q , r ) s 0. P ( q, r ) = q 2 - qr + r'\ P^tq.r) = 2q - 3q r - 3qr + 2r . etc. ) 
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is constructed according to some general rule (the "rule" here is the required 

q, r-symmetry) , there comes a point in the construction where 1 be set of all 

available independent such forms is exhausted and any form constructed 

thereinafter contains absolutely no new information y being merely a 

modification/repetition of what lias already preceded. What this means for the 

case at hand is that the most general sixth-order form must be a linear 

3 2 

combination of the products. F^(q,r) and fq,r), while a seventh-order such 
form is { F *( q, r )F^(q, r)} . A simple appropriate algebraic expression for 
Q(q,r\t) which for q s 0 fits data reported by Stewart^ is thus 

AF^fq, r ) 

Q(q,r\t) « cr 3 (t) - 

{1 + BF (q.r)P 
l ^ 

~3 —2 

where A = /3A (t), B = ^A (t), and f3,y ~ dimensionless constants; 



The specific (3 and y for the data sets highlighted by Stewart (loc. cit.) are: 
Curve 1 (cf. Figure 7); 13 = -0.74, y = 7; 

Curve 2 (cf. Figure 8): 13 - - 1.5 , j = 12: 

Curve 3 (cf. Figure 9); 0 = -3.71, y = 22: 

The value of C^ which most closely approximates all three cases is C - 
- 0 . 08 tt. 


6. Conclusions and Concluding Remarks 

Since integral scale provides a measure of turbulence randomness 
realistic estimates of this scale are essential to ensure related flight 
simulation fidelity. The role played by anisotropy in a wind shear encounter 
is quantified in the correlation length, L = A( t j + AA( t ) . As an aircraft 
penetrates the headwind portion of Shear Flow I. AA > 0 since M and a are 


15 



negative in this phase of flight. The turbulence sensed there by the aircraft 


is therefore less random than that predicted by isotropic theory. The 
tailwind phase of such an encounter provides AA < 0 — a more random turb ul ence 
in comparison with isotropy. A similar analysis of Shear-Flow II yields 
negative AA on both sides of the microburst, thus making the attendant 
turbulence more random than the isotropic case in all phases of the encounter. 

The significance of anisotropy in microburst turbulence models has herein been 
considered from the perspective of power spectral density. The energy scaling 
factor clearly defines the inadequacy of using isotropic theory only to 
estimate the turbulent energy in a wind shear. The psd representat ion of' the 
turbulent downwash is a step towards enhanced simulation of low-level 
turbulence, with the relative importance of the characteristics considered 
above becoming clear only upon eventual flight-simulator implementation. 

Acknowledgements 

This research was supported by the NASA-Langley Research Center through the 
ASEE and LARSS programs. The authors wish to thank Drs. R. L. Bowles and B.T. 
McKissick (both of NASA-LaRC) and Mr. D. Staublin and Mr. A. Rosenberger (both 
of Michigan Technological University) for their collective participation in 
this work. 


16 


References 


1 ... T T "Spearhead Echo and Dovmburst Near the Approach End of a 

F u,nta, T T ■ , Spear head MesometeoroloKicai Research Project, 

John F. Kennedy Airport , batei_LL!^ siiSi * — - 

SMRP Research Paper 137, 51 PP- , 1976. 

‘'Taylor, 0 ... "Statisticai Theory of Tur b ulen=e, '' EC****^ « U» 

Royal Society ;, A, Vol. 151, 1935, PP- 

o . . . . T "On the Statistical Theory of Isotropic 

3 von KcLrm&n, T. and Ho wart h, L. , . v , 164 1937, pp. 192-215. 

Turbulence," Pr^eedin^ of the Royal Society, A, Vol. 164, 

^Robert son H.P., "Invariant Theory of Isotropic Turbulence^ Proceedings 

pf tJ^^ Vo1 - 36 - 1940 ’ pp - 209 ' 233 - 

Batchelor. G. K. , "Theory of Axisymmetric Turbulence, " Proceed ings o_f 

the RoyaT Society;, A, Vol. 186, 1946, pp. 480-o02. 

®Trevj no, G. , "Turbulence for Flight Simulation," .iournal of Aircraft, 

Vol. 23. April 1986, pp- 348-349. 


7„ 


Treviho. G. . "Airplane Flight Through Hind-Shear Turbulence." Journal o£ 
Aircraft , Vol. 23, Sept. 1986, pp. 733 735. 


a Qt rnr f ur ^ of Nonstationary Atmospheric 

^Trevino, G. . T> me-' n»n .ant t Stn ucture 827 - 828 . 

Turbulence," Journal of Aircraf l, Vo ■ 2l -- p 

Batchelor, G. K. , Theory of Homogeneous Turbulenc e, Cambridge University 
p re ss, Cambridge, U. K. , 1967. 

IH Trevl«o. G. , "Turbulence Structure in Microburst Phenomena," Journal, 
of Aircraft , Vol. 24. April 1987, pp. 283 285. 

11 , _ . . u u "On the Decay of a Normally Distributed 

1 Proudman, I. and Reid, • > ,, Phi | osophical Transactions of the 

and Homogeneous Turbulent : Verity ^ Field 

Royal Society of London, Vol. 247. A. 92b, pp 

12 Tr , vtno c„. "invariant Structure In Nonstationary Behavior", 

forthcoming in Journal sL SoMEd and Vibrat ion 


Hilbert. D. , "Ober die Theorie der Algebraischen Formen. " Mathematisc he 
Annalen, Vol- 36 . 1890 , pp. 473 - 534 . 

1 ^Stewart . R.H.. "Triple Velocity Correlations Jn ^tropic Turhnlenc<r" 

Proceedings. .of the Bunbridgn Pin losophi . c ai Scoie_o. 

I 46-157. 


13. 


17 


Appendix: The Case for Self-Preservation 


Nonstationary behavior whose time-evolution can be effected by two 

time-dependent scaling effects is important not only because of its relatively 

simple features, but also because of its immediate practical application. For 

nonstationary processes the two-time correlation function, QU , t ) = 

<w(t )w(t )>, can be considered as a two degree-of-f reedom system in which I 
1 2 

= (t + t 1/2 and t = t - t play the role of uncoupled principal 

12 2 1 

coordinates; indeed in the ( t . x ) -coordinate reference frame the two 
degree-of-f reedom system at hand is instead duly converted into one whose 
physical characteristics accordingly vary with time, i.e. Q(t *t ) a Q(t,x) a 
C ( t , r/V ) - C(t,r), where V * rx. And even though time-varying systems can be 
regarded as filters which "select” eigenfunctions in the pure sense rather 
than just simple sinusoids, the natural tendency in their analysis is to treat 
them as systems whose properties do not change "appreciably" with repeat to 
time. The underlying motivation for this approach is that classical modal 
analysis is not explicitly restricted to systems which are time-invariant 
only; it is only necessary that the differential equation governing the 
system, as well as the boundary conditions, be separable into a function of 
time alone and in this case a function of non-dimensional ized length alone, 
i.e. that temporal changes in the system be in some sense "uniform". The 
suggested separation for dynamical systems is thus the essence of 

self-preservation for nonstationary behavior in random processes. The 

mathematics of same proceeds as follows. Consider the most general case of 
time-varying intensity and integral scale, viz. C( t . r |=<r( t-x/2 )<r( t -x/2 )Rf t , t I . 
r = r/V. if' cH * * - ) is expanded in Taylor's series about x = 0. the result for 
r: ( t , X ) i s 

m.T ) - /( U { 1 - t 2 /4TO R( t . t | 


10 


where only linear terms in t have been retained for <r( t ± r/2 ) , and T 

is the time-dependent scale over which <r(t) may be assumed to vary 

T=0 

Linearly with t, viz. 



<r( t ± x/2 ) - cr( t. ) { 1 ± x/2T} . 
Since R( • • * ) scales with Alt), i.e. R(t,x) ~ Hlx/A) 
T, the modified expression. 


and <r( * * * ) scales with 


C( t , x) » <r 2 ( till - (x/2A) 2 ( A/T)"}Rl - - - ) , 
allows for analysis of the coupling between the two scales, T and A. Note 
that as far as R( * • • ) is concerned, R( * • • ) - 0 for x -> A. so that the ma jor 
contribution to C(t,x) from Rl * - ■ ) itself is from values of x < A. Therefore, 
if T » A, it follows that Clt.x) - <r 2 ( t )R( * • * ) , i-e. the x-dependence of 
C(t,x) is determined by R( ■ * * ) only and is unaffected by the x dependence of 
{<r(t - x/2)«r(t + x/2)> 


*A Alt) in this expression considered as a "time length scale 
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Ilgure 1: 

F i gure 2 : 

F i gure 3 : 

Figure 4; 

Figure 5: 
Figure 6: 

Figure 7: 

Figure 8 

Figure 9 


Figure Captions 

Schematic of a typical microburst encounter by an aircraft 
attempting to land. 

Contributions of isotropy and anisotropy to the downwash 
corre lat ion. 

Time-evolution of a self-similar or self-preserving correlation 
function; note the explicit two scaling-effects feature. 

Energy scaling factor quantifying the effects of anisotropy, i.e. 

A A/ A. 

Sheared microburst mean wind consistent with no slip 

Sheared microburst mean wind for continuously varying 

headwi nd/t ai 1 wi nd . 

Triple-velocity curve. *<"•>. obtained by Stewart, broken line is 
theoretical curve. 

Triple-velocity curve. obtained by Stewart; broken line is 

theoretical curve. 

Triple-velocity curve. k<->, obtained by Stewart; broken line is 
theoret. ical curve 
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experimental curve 2 
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experimental curve 3 
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